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We present selected results of our program to determine the masses, gauge cou- 
plings, and Yukawa couplings of the minimal supersymmetric model in a full one- 
loop calculation. We focus on the precise prediction of the strong coupling Oc s (Mz) 
in the context of supersymmetric unification. We discuss the importance of includ- 
ing the finite corrections and demonstrate that the leading-logarithmic approxima- 
tion can significantly underestimate Ct s (Mz) when some superpartner masses are 
light. We show that if GUT thresholds are ignored, and the superpartner masses 
are less than about 500 GeV, the prediction for OL s {Mz) is quite large. We impose 
constraints from nucleon decay experiments and find that minimal SU(5) GUT 
J> \ threshold corrections increase Qt s (Mz) over most of the parameter space. We also 

consider the missing-doublet SU(5) model and find that it predicts preferred values 
\ for the strong coupling, even for a very light superpartner spectrum. We briefly 

■ discuss predictions for the bottom-quark mass in the small tan p region. 

o 
in 

^ ■ 1. Introduction 

Q_i! The exact one-loop corrections to the masses, gauge couplings and Yukawa cou- 

plings of the minimal supersymmetric model are described in Ref. [Q. These correc- 
tions are essential ingredients for accurate tests of grand unification. They allow one 
to extract the underlying TJR parameters from a given set of measured observables. 
The ~DR parameters can then be run up to a high scale to explore the consequences 
^ ■ of different unification hypotheses. 

Alternatively, the radiative corrections can be used to translate various limits into 
excluded regions of the TJR parameter space. This is illustrated in Fig. 1, where 
we show the excluded region of the Mo, M1/2 parameter space at the tree and one- 
loop levels, from current experimental constraints.* Here, "tree- level" means that 
the superpartner masses are determined from the TJR parameters evaluated at the 
scale Mz- For the Higgs mass, the tree- level curve corresponds to the one-loop mass 
neglecting the gauge/Higgs/gaugino/Higgsino contributions. 

The study of gauge coupling unification has been carried out by many groups 
beginning more than 20 years ago. There has been a resurgence during the last 
five years. The analyses have become increasingly more refined. The most recent 
analyses use the full set of two-loop RGE's to predict the strong coupling constant as 
a function of the electroweak input parameters. They pay proper attention to MS - 
DR differences and treat the weak thresholds in different levels of detail. 



"Mq is the universal scalar mass, My 2 is the universal gaugino mass, and Aq is the universal A-term. 
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Figure 1: Excluded region (shaded) of the M , Mi/ 2 plane, for tan j3 = 2, m t = 170 
GeV, and A$ = 0. All masses are evaluated at one-loop. The symbols indicate which 
experimental constraint is relevant: \ + =^ m x + ^ ^ GeV; g m g > 125 GeV; 
v => rriy > 42 GeV; h => > 58 GeV; x° direct searches for neutralinos at 
LEP. The dashed line shows the boundary of the excluded region at tree-level. 



Over time, the predicted value of the strong coupling constant has increased 
markedly, in part due to the refinements mentioned above, but more so from the 
fact that the standard-model weak mixing angle, as determined by a global fit to the 
data, has been steadily decreasing. (This is correlated with the increasing best-fit 
value for the top-quark mass.) 

In this talk we take a closer look at supersymmetric unification. We treat the 
super symmetric threshold corrections in a complete one- loop analysis.^ Our work 
stands in contrast to most previous studies, which are based on the "leading loga- 
rithm approximation." This approximation involves taking the standard-model value 
of sin 2 6w and adding the logarithmic parts of the SUSY threshold corrections. The 
approximation works well if all of the SUSY particle masses are much greater than 
Mz, in which case the decoupling theorem implies that the finite effects of the SUSY 
particles are negligible for all low-energy observables. 

However, in realistic models it is not unusual for the supersymmetric spectrum to 
contain light particles of order the Z-mass. In this case one cannot use the standard- 
model value of sin 2 9w as an input into a precision analysis. This is because the quoted 
value of sin 2 9w is the result of a fit to the data, assuming that the standard model is 
correct. The experimental analyses do not include the finite SUSY corrections, which 
are different for each observable. Therefore in our analysis, we use a single set of 

tSee Ref. [§ for a similar treatment of finite corrections to sin 2 9w ■ 
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Figure 2: Contours of a s (Mz) in the M , M 1 / 2 plane, with tan/? = 2, m t = 170 
GeV, and Aq = 0. The contours of squark mass 500 and 1000 GeV are shown by 
dot-dashed lines. The x's in the upper left hand corners indicate regions in which 
the tau sneutrino is the LSP. The shaded regions are excluded by particle searches. 



inputs in our calculation of sin 2 6 W , namely, ocem, Gp, M z , m t , and the parameters 
that describe the supersymmetric model. 

We note that a careful evaluation of the weak mixing angle is important for 
determining a precise prediction for a s (M z ). Using the one-loop RGE's and the 
condition of coupling unification, we find that the three gauge couplings satisfy 

fo-A fa-fa m 
g 2 M fifr) $(a0 " ' u 

where 0i are the three beta functions, dgi/dt = 0igf/167r 2 . These relations imply that 

a s (M z ) ~ a Lb P ■ {Z) 

Hence, an error in the determination of s 2 of 1% leads to an error in a s [Mz) of 7.5%. 

In this talk we use the full set of one-loop radiative corrections to evaluate the 
UE gauge and Yukawa couplings. The ~DR couplings serve as the boundary condi- 
tions for the two-loop gauge and Yukawa coupling renormalization group equations 
(RGE's), which determine the couplings at very high scales. In what follows we use 
the full one-loop corrections at both the weak and GUT scales to determine the re- 
gions of supersymmetric parameter space that permit gauge and Yukawa coupling 
unification. 
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Figure 3: Predictions for a s {Mz) using three different values for the QED coupling 
a(Mz)- In each case we show the ±1-<t bands, for tan/? = 2, m t = 170 GeV, 
Aq = 0, and \i > 0. 



2. Calculation of s 2 

Given the inputs a EM = 1/137.036, G F = 1.16639 x 1(T 5 GeV -2 , and M z = 
91.187 GeV, as well as m t and the parameters of the supersymmetric model, we 
determine (c 2 = 1 — s 2 ) 

OtEM j -2 -2 17 a EM 

a = tz r-7T and c s — 



(1-Ad) ^ V2G F M Z {1- Af) ' 

Here Ad contains logarithms of the masses of the charged particles, and 

. . flz(M 2 ) IL w (0) 

Af = Ad ; f9 z; + ™\ ' + vertex + box . 

M 2 M 2 , 

The vertex and box diagram contributions are the so-called "non-universal" or "non- 
oblique" corrections, and the remaining corrections involve the real and transverse 
parts of the ~DR gauge boson self-energies. The vertex and box corrections vanish 
in the leading logarithm approximation; the correction Ad contains only logarithms; 
and the W and Z self-energies contain both logarithmic and finite corrections. 

In our calculation of Af we include the dominant two-loop corrections (given in 
Ref. Jl), which leads to a very precise determination of s 2 . Following Ref. 0, we 
estimate the theoretical uncertainty in s 2 to be about 1 part in 10 4 , while the exper- 
imental uncertainty (due to the uncertainty in the determining the electromagnetic 
coupling at the Z-scale) is 2.6 parts in 10 4 . Having determined d and s 2 precisely, 
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Figure 4: Comparison of a s {Mz) in the leading logarithm approximation (LLA) 
versus the full one-loop calculation, for tan f3 = 2, mt = 170 GeV, Aq = 0, and fi > 0. 
The dashed line shows the result if the non-universal corrections are neglected. 



we are in position to fix the boundary conditions for the two- loop RGE's ||, 




and accurately investigate gauge coupling unification. 

In the following we assume that the SUSY masses unify at the scale Mqut (which 
is defined as the scale where gi and §2 meet). Therefore the supersymmetric model 
is parametrized by a universal gaugino mass M1/2, a universal scalar mass Mq, and 
a universal trilinear scalar coupling Aq. The ratio of ~DR vacuum expectation values 
evaluated at the scale Mz is denoted tan/3. We require the parameters to be such 
that electroweak symmetry is spontaneously broken. These conditions determine the 
value of fi 2 , the supersymmetric Higgs mass parameter, and tiia, the CP-odd neutral 
Higgs boson mass, once we specify the sign of /i. 



3. Prediction for a s (M z ) 

As a reference point, we show in Fig. 2 contours of a s (Mz) in the M , M1/2 
plane, with no GUT thresholds, tan/5 = 2, m t = 170 GeV, and A Q =0. We find 
a s (M z ) is large compared to the PDG value |J a s (M z ) = 0.117 ± 0.005, especially 
near M , M 1/2 = 100 GeV. 

The experimental uncertainty in the determination of a s (Mz) is primarily due 
to the uncertainty in determining the electromagnetic coupling at the Z-scale. We 
use the value recently determined by Eidelman and Jegerlehner 0. Martin and 
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(b) M 1/2 = 80 GeV 

I I I I I I I I I I I I I I I I I I I I I I I I I 

50 100 200 500 1000 2000 
M (GeV) 



Figure 5: The smallest triplet Higgs mass allowed by nucleon decay constraints in 
minimal SU(5), with tan/3 = 2, m t =170 GeV, A = 0, and fi > 0. The GUT scale 
is also shown. 



Zeppenfeld |§ and Swartz |§ have also performed analyses to determine a{Mz). We 
show in Fig. 3 the differences in the determination of a s (Mz) using the various values 
oia(Mz). 

As stated in the introduction, the finite corrections can be significant when some 
of the superpartners have masses of order Mz- This is illustrated in Fig. 4, where we 
compare the value of a s (Mz) in the leading logarithm approximation (LLA) with the 
value obtained in the full calculation. In Fig. 4(a) the full and LLA curves converge 
for large Mi/ 2 because the SUSY particles decouple. In Fig. 4(b) the full and LLA 
curves do not converge as M becomes large. This is because Mi/ 2 = 60 GeV, so the 
gauginos remain light for arbitrarily large Mq. 

To summarize our results for a s [Mz) m the absence of GUT threshold corrections, 
we find a s (Mz) > 0.126 for squark masses less than 1 TeV, with m t = 170 GeV. For 
a SUSY spectrum of 500 GeV or less, we have a s {Mz) > 0.130. 

If we require smaller values of ot s {Mz) and a light supersymmetric spectrum, a 
GUT threshold correction is clearly needed. We can parametrize the GUT threshold 
correction by e g , where 

<7 3 (M GUT ) = <?gut(M GU t) (1 + e g ) , 

and <?gut = <?i(Mgut) = ^(Mgut)- A smaller value of a s (Mz) requires e g < 0. In 
what follows we examine the value of e g in two SU(5) GUT models. 



In the minimal SU(5) model fll0" |, the gauge coupling threshold correction e' is 
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Figure 6: Contours of the smallest possible a s {Mz) consistent with nucleon decay 
in minimal SU(5), with tan/3 = 2, m% = 170 GeV, and Aq = 0. The dot-dashed 
lines indicate contours of 500 and 1000 GeV squark masses. The x 's indicate the 
region where the LSP is charged, and the lower shaded regions are excluded by 
particle searches. 



given by |TT 



40tt 2 ° \M GVT J 

where M# 3 is the mass of the color-triplet Higgs particle that mediates nucleon decay. 
From this expression, we see that e' g < whenever M# 3 < Mqut- However, M# 3 
is bounded from below by proton decay experiments. The M# 3 mass limit is of the 
form Wl 



\l+v tK \ 
Mh 3 > M ' . y — l f{w, d, u, e 
smzp 

where M. is a nuclear matrix element, y tK parametrizes the amount of third generation 
mixing, and / is a function of the wino, squark and slepton masses. 

In Fig. 5 we show the minimum value for Mh 3 from the nucleon decay constraint, 
for the conservative choices M. = 0.003 GeV 3 and |1 + y tK \ = 0.4. We see that 
Mf 3 in > M GUT unless M > 500 GeV and M 1/2 < M . Thus, in most of the parameter 
space, e' g > 0. 

In minimal SU(5), a s (Mz) is typically even larger than it was without any GUT 
thresholds, as illustrated in Fig. 6. The only exception occurs in the region M 3> 
Mi/2, where the proton decay amplitude is suppressed. In this region, with 1 TeV 
squark masses and m t = 170 GeV, we find a s (Mz) as small as 0.123. In fact, as long 
as rriq < 1 TeV, a s (Mz) < 0.126 can only be obtained in the region M ~ 1 TeV. For 



7 



Missing doublet SU(5) 



/x > fi < 




100 200 500 1000 2000 100 200 500 1000 2000 



M (GeV) M (GeV) 

Figure 7: Contours of the minimum possible a s (Mz) consistent with nucleon decay 
in the missing-doublet model, with tan/3 = 2, mj=170 GeV, and Aq = 0. The 
shading and dashed lines are as in Fig. 6. 



example, if M < 500 GeV, a s {M z ) > 0.127. 

The missing-doublet model is an alternative SU(5) theory in which the heavy 
color-triplet Higgs particles are split naturally from the light Higgs doublets [y|. In 
this model the GUT gauge threshold correction is given by Jl4| 



<' = IS? { ,og S) - T log 5 + 15 ,og 2 } = " 4% ' (4) 

Thus, for fixed M# 3 , the missing-doublet model has the same threshold correction as 
the minimal SU(5) model, minus 4%. In eq. (|J), Mjj is the effective mass that enters 
into the proton decay amplitude, so the bounds on Mh 3 in the minimal SU(5) model 
also apply to Mjj in the missing-doublet model. 

The large negative correction in eq. (ffl) is due to the mass splitting in the 75 
representation, and gives rise to much smaller values for a a (Mz). This is illustrated in 
Fig. 7, where we show contours of a s (M z ) in the M , Mi/ 2 plane, with = M^ n . 
The values of a s (Mz) are somewhat low, but one can easily obtain larger values, for 
example, by increasing Mfj . 

In Fig. 8 we illustrate the full range of a s (Mz) values in the minimal SU(5) and 
missing- doublet models. We show the allowed range of e g in the two models (setting 
M^ 3 ax = 10 19 GeV), together with the range of e g which yields a s {M z ) = 0.117±0.01. 
The missing doublet model is mostly contained within the preferred region of a s (Mz), 
while minimal SU(5) is almost entirely outside. 
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Figure 8: The light shaded regions indicate the allowed values of the gauge coupling 
threshold correction e g in the minimal and missing-doublet SU(5) models. The dark 
shaded region indicates the range of e g necessary to obtain a s (Mz) = 0.117 ± 0.01. 
For tan/? = 2, m t = 170 GeV, A = 0, and p > 0. (From Ref. fl|.) 



4. Yukawa unification 

In the final part of this talk, we investigate the possibility of bottom-tau Yukawa 
coupling unification. Our procedure is as follows. We start with the experimental 
value for the r pole mass m T = 1.777 GeV ||. We convert it to the DR running mass 
and evolve it up to the Z-scale, where we apply the SUSY threshold corrections. We 
compute the ~DR vev from the Z-boson mass, v 2 = 4(M| + fl z (Mz))/(g' 2 + g 2 ), and 
use it to determine the UR tau Yukawa coupling 

- V2rh T 

A r — n ■ 

v cos p 

We then solve the RGE's to find A t (M G ut), and set the bottom Yukawa coupling to 

VM GUT ) = A T (M GUT ) (1 + e b ) , 

where e& parametrizes the GUT threshold correction. Once we have A&(Mgut) we 
run everything back to the weak scale and self-consistently determine the pole mass 
for the bottom quark. 

Let us first examine the prediction for the bottom-quark pole mass with e& = 0. 
Generally, the large value of the strong coupling increases the bottom mass so much 
that the prediction typically falls outside the region determined by experiment ||, 
which we take to be 4.7 GeV < m h < 5.2 GeV. Because Yukawa couplings enter the 
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Figure 9: The bottom-quark mass and a s {Mz) vs. mt for the case of no GUT-scale 
thresholds, for various values of tan /?, with Aq = 0, \i > 0, and Af (Mqut) = 3. The 
right (solid) leg in each pair of lines corresponds to varying from 60 to 1000 

GeV, with Mq fixed at 60 GeV. The left (dashed) leg corresponds to Mq varying from 
60 to 1000 GeV, with M]/ 2 = 100 GeV. On the solid lines the circles mark, from 
top to bottom, M 1/2 = 60, 100, 200, 400, and 1000 GeV, and on the dashed lines 
the circles mark M = 60, 200, 400, and 1000 GeV. Note that the lowest point on 
each left leg and the second-to-lowest point on each right leg corresponds to rriq ~ 1 
TeV. The horizontal dashed lines indicate mj, = 5.2 GeV and a s (Mz) = 0.127. The 
x's mark points with one-loop Higgs mass mh < 60 GeV. (Prom Ref. fi~5| .) 



Yukawa RGE's with the opposite sign from gauge couplings, large Yukawa couplings 
help reduce the large b mass. In particular, in the very small tan/? region, the top 
Yukawa coupling becomes large. In this infrared fixed point region the bottom mass 
can be less than 5.2 GeV. 

We show in Fig. 9 the prediction for m& and a s (Mz), obtained by setting the top 
Yukawa coupling to be A 4 (M G ut) = 3. We see that even with this large top Yukawa 
coupling, which is on the verge of being non-perturbative, m& is larger than 5.2 GeV 
unless M or Mi/ 2 is greater than about 1 TeV. One needs a GUT threshold correction 
to reduce m b with a SUSY mass scale below 1 TeV. 

The similarity between the curves for a s [Mz) and m& in Fig. 9 illustrates the 
strong correlation between a s (Mz) and mj,. In fact, m& is far more sensitive to the 
gauge-coupling GUT threshold correction than that of the bottom-quark Yukawa 
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Figure 10: The bottom-quark mass vs. the top-quark mass for fixed values of 
a s (Mz) and various tan/3. The solid lines correspond to A^Mqut) = 3 while the 
dashed lines correspond to A^Mqut) = 2. The upper line in each pair corresponds 
to a light super symmetric spectrum with Mq = Mi/ 2 = 80 GeV. The lower line in 
each pair corresponds to a heavy spectrum, Mo = 1000 GeV, M 1 / 2 = 500 GeV. The 
dotted lines delineate the preferred region 4.7 < m\> < 5.2 GeV. For ^0 = and 
H > 0. 



coupling. Numerically, we find 

Sin 

— ~ 0.8e b + 8s Q . 
m b 

Hence, if we consider a GUT model where e g is sufficiently negative, the central value 
of mb can De obtained with e h = 0. This is shown in Fig. 10, where, for fixed e h = 0, 
we show the predicted value of nib, assuming various values of e g that yield particular 
values of a s (M z ). We show the results for A t (M GUT ) = 2 and 3, and for a small and 
large supersymmetric mass scale. The figure shows that as long CIS Eg IS such that 
a s (M z ) — 0.12, an acceptable value of m b is predicted for e b = 0, independent of the 
top-quark mass and the supersymmetric mass scale. 

5. Conclusion 

In this talk we have presented results from a complete calculation of the one-loop 
corrections to the masses, gauge, and Yukawa couplings in the MSSM. We have seen 
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that such a calculation allows us to reliably investigate various unified models to see 
whether they are compatible with current experimental data. 

In particular, we found that the finite SUSY corrections, which are neglected in 
the leading logarithm approximation, can substantially increase the prediction for 
a s (Mz) when some of the SUSY partner masses are lighter than or of order Mz- In 
the minimal SU(5) model, we found that a s (Mz) ^0.14 in the small Msusy region, 
M ~ M1/2 — 100 GeV. We also found a s (M z ) > 0.123 in the region where the squark 
masses are below 1 TeV. In contrast, we showed that the missing-doublet SU(5) 
model can accommodate much smaller values of a s (Mz), such as a s (Mz) — 0.113 for 
M ~ M 1/2 ~ 100 GeV. 

This work was supported by the U.S. National Science Foundation under grant 
NSF-PHY-9404057. 
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